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We study the Kondo effect in a quantum dot coupled to two noncollinear ferromagnetic leads. 
First, we study the spin splitting 5e = ej — e-f of an energy level in the quantum dot by 
tunnel couplings to the ferromagnetic leads, using the poor man's scaling method. The spin 
splitting takes place in an intermediate direction between magnetic moments in the two leads. 
8e oc p\/cos 2 (6>/2) + v 2 sm 2 (6>/2), where p is the spin polarization in the leads, 9 is the angle be- 
tween the magnetic moments, and v is an asymmetric factor of tunnel barriers (— 1 < v < 1). Hence 
the spin splitting is always maximal in the parallel alignment of two ferromagnets (9 = 0) and 
minimal in the antiparallel alignment (8 — n). Second, we calculate the Kondo temperature Tk. 
The scaling calculation yields an analytical expression of Tk as a function of 6 and p, Tk(9,p), when 
Se <C Tk. Tk(9,p) is a decreasing function with respect to p\/cos 2 (6/2) + v 2 sin 2 (#/2). When Se is 
relevant, we evaluate Tj<.(Se,8,p) using the slave-boson mean-field theory. The Kondo resonance is 
split into two by finite Se, which results in the spin accumulation in the quantum dot and suppression 
of the Kondo effect. 

PACS numbers: 75.20.Hr, 72.15.Qm, 72.25.-b, 73.23.Hk 
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I. INTRODUCTION 

The Kondo effect, a typical many-body effect in metals 
with dilute magnetic impurities, has given exciting top- 
ics for decades i The Kondo effect gives rise to the for- 
mation of local singlet states between itinerant electrons 
and localized spins. A characteristic energy scale, Kondo 
temperature Tk, corresponds to the binding energy of 
the singlet states. At temperature T < Tk, the fluctu- 
ation of localized spins is significantly enhanced by the 
spin-flip scattering of conduction electrons at the Fermi 
level. Simultaneously, the scattering of the conduction 
electrons is magnified, and in consequence the resistance 
increases in bulk metals. 

Recent progress in fine processing technologies has en- 
abled us to make artificial atoms on semiconductors, 
quantum dots (QDs). In QDs, discrete energy levels are 
occupied by a fixed number of electrons by the Coulomb 
blockade. With an odd number of electrons, the QDs 
behave like magnetic impurities. The Kondo effect has 
been observed when they are coupled to external leads. 2,3 
In this case, the Kondo effect enhances the conductance 
through the QDs since the conduction electrons are res- 
onantly transported through the Kondo singlet state. 

Lately, spin-dependent transport has been of great in- 
terest from a viewpoint of spin-based electronics, so- 
called spintronics^ A notable example is the tunnel mag- 
netresistance in a spin- valve geometry, in which two ferro- 
magnetic leads are separated by a non-magnetic insulat- 
ing layer £ The tunnel current is proportional to pl-\-p, in 
the parallel (P) alignment of two ferromagnets, where pf 
and pi are the density of states for majority and minor- 
ity spins, respectively. It is proportional to 2p-\p\, in the 
antiparallel (AP) alignment, which is smaller than that 
in the P alignment for p^ > p±. When the P alignment is 



changed to AP alignment by applying a magnetic field, 
an extremely large magnetresistance is observed.^' 8 - The 
spin-valve effect in a QD sandwiched by two ferromag- 
netic leads has also been investigated^^ 

In this article, we pay attention to the 
Kondo effect in the spin-valve geometry includ- 
ing a QD. This Kondo effect has been studied 
theoreticallyii^^ 3 -^^^^^^^^ with motives 
elucidating how the ferromagnetism of the leads influ- 
ences the fluctuation of a localized spin. One of the 
theoretical problems was whether the energy levels in 
the QDs are spin-spli t 13 ' 15 or no t 19 ' 20 in the absence of 
magnetic field. The numerical renormalization group 
technique, which is one of the most reliable methods 
to treat the Kondo physics, concluded the existence of 
the spin splitting, Se = — ef , 16 ' 17 The spin splitting 
arises from the spin-dependent charge fluctuations by 
the tunnel coupling to the ferromagnets and weakens 
the Kondo effect. 18 The spin splitting has been observed 
in experiments as the suppression of zero-bias anomaly 
of the differential conductance, which is a hallmark of 
the Kondo effectj 22 ' 23 The splitting is the largest in 
the P alignment and smallest in the AP alignment, in 
accordance with the theoretical results. 

Until now, the Kondo effect in QD spin-valves 
has been elucidated for collinear cases of two fer- 
romagnetic leads, P or AP alignments, by numer- 
ous papers . 12 d 3 i 15 i 16 ] 17 i 18 i 19 i 20 However, there exist 
very scarce theoretical investigations for noncollinear 
alignments ! 11 ' 14 ' 21 The purpose of this article is to sys- 
tematically study the Kondo effect in QDs coupled to 
noncollinear ferromagnetic leads. The spin splitting and 
Kondo temperature Tk are evaluated for arbitrary align- 
ments of the ferromagnetic leads. 

In the first part of this article, we generalize the dis- 
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cussion based on the scaling approach adopted for P and 
AP alignments in Ref. [l3|. For the Kondo physics, all the 
energies from Tk to the upper cutoff (band width Do) 
have to be taken into account properly. We use the poor 
man's scaling procedure^ in two stages^ First, we re- 
duce the energy scale D until the charge fluctuation is 
quenched [D = D\ given by Eq. (|21[) ] and renormalizc 
the energy levels in a QD (first stage scaling). Then 
we obtain the spin splitting 5e, which is directed in an 
intermediate direction between the magnetic moments in 



the two leads. Se oc py cos 2 (9/2) + v 2 sin 2 (9/2) in the 
absence of magnetic field, where p is the spin polariza- 
tion in the leads, 9 is the angle between the magnetic 
moments, and v is an asymmetric factor of tunnel barri- 
ers (-1 < v < I). 

We proceed to the second stage of the scaling for D < 
Di, where a localized spin fluctuates in the QD with 
fixed spin splitting Se. The exchange couplings for the 
localized spin increase with decreasing D until the scaling 
stops at D = Se or reaches a fixed point of the strong 
coupling limit corresponding to the Kondo effect (D = 
Tk). When Se <C Tk, the Kondo effect takes place. We 
obtain an analytical expression of Kondo temperature Tk 
as a function of 9 and p. Tk(#, p) is a decreasing function 



with respect to p^J cos 2 (9/2) + v 2 sin 2 (9/2). When Se > 
T K , no Kondo effect is expected. 

When Se is comparable to Tk, the scaling method does 
not work. In the second part of this article, we evalu- 
ate T K (Se, 9,p) using the slave-boson mean-field (SBMF) 
theory^ The SBMF theory describes the Kondo reso- 
nant state on the assumption of its presence and Fermi 
liquid behavior. The theory is exact for the Kondo effect 
in (N — oo)-fold degenerate Anderson model with infi- 
nite U in the case of non-magnetic leads. For N = 2 (spin 
S = f/2), it is widely used for the semi-quantitative esti- 
mation of the Kondo temperature, the Kondo resonance 
at the Fermi level, and resonant transmission. However, 
we have to take a special care to apply the SBMF to the 
present case with magnetic leads. In the slave-boson for- 
malism, the empty state in the QD is expressed by an 
auxiliary boson field b and singly-occupied states with 
spin a =Til are by pseudo-fermion operators f a . The 
SBMF theory describes the Kondo effect by the conden- 
sation of the boson field, replacing the boson operator 
b to be a c number. We observe the formation of the 
Kondo resonant state at the Fermi level with the width 
of Tk- In this mean-field theory, the spin fluctuation (de- 
scribed by f a ) is considered to construct the Kondo res- 
onance on one hand, the charge fluctuation is neglected 
(b is replaced by a constant) on the other hand. Indeed 
the spin splitting Se cannot be evaluated by the SBMF 
theory . 12 ! 20 To discuss the Kondo effect induced by the 
spin fluctuation, we apply the SBMF theory to the renor- 
malized Hamiltonian with the band width D\ after the 
spin splitting Se by the charge fluctuation has been de- 
termined in the first stage scaling. 

When Se = 0, the SBMF theory yields the for- 



mation of the Kondo resonant state at the Fermi 
level. The Kondo temperature, which is defined by 
the width of the resonance, decreases with increasing 

p\j cos 2 (9/2) + v 2 sin 2 (#/2), in semi-quantitative agree- 
ment with that obtained by the scaling method. This 
indicates an applicability of the SBMF theory to the 
present problem. A finite Se splits the Kondo resonance 
into two and suppresses the Kondo effect. We evaluate 
Tk(5e, 9,p) by the width of the split resonances. We also 
calculate the conductance G at temperature T = 0, which 
is determined by the resonant tunneling through the 
Kondo resonant state. G could show a non-monotonic 
behavior as a function of Se, in contrast to Tk- This 
is ascribable to an interference effect between two spin 
components, analogous to the Fano resonance. 

This article is organized as follows. In the next section 
(Sec. II), our model is presented. In Sec. Ill, the spin 
splitting Se of the QD level and the Kondo temperature 
Tk are calculated using the poor man's scaling method. 
In Sec. IV, we examine Tk and conductance by the SBMF 
theory. The conclusions are given in Sec. V. 



II. MODEL 

We consider a single-level QD coupled to two ferro- 
magnetic leads as shown in Fig. [TJ The ferromagnets are 
identical and aligned in a noncollinear way: The mag- 
netic moment in the lead L (R) is tilted by 9/2 (—9/2) 
from the z axis in the z-x plane (0 < 9 < tt). This model 
is described by the Anderson Hamiltonian, 

H = ^2 e k<7C Lk la c L k,i<T + ^2 ek*c Rk rCT Cij fc , rCT 
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+ X^cL , ff di«, + V R J Rk d ra + h.c), (1) 
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where CLk.io- (cRk.ra) is a fermion operator of an electron 
with wavenumber k and spin la (ra) in lead L (R), and 
d za is that with spin za in the QD level. We denote 
^T/^l( r T/ r l) f° r majority /minority spin in lead 
L (R) and z\ /z \ for spin-up/down in the z direction 
in the QD. Fermion operators in the QD, dia, d ra , are 
related to d za by 

di\ 
d U 

d„ 



cos 3 4 T - 
- sin 1 4 



cos 4gL 



smja zl , 



rt 



sin ■ 



dri = sin jd z i 



cosjd zl . 



(2) 
(3) 



The density of states in the ferromagnetic leads is con- 
stant, and pi for majority and minority spins, respec- 
tively (pf > pi), in the band of —Do < uj < Do. The 
spin polarization of the leads is given by 



PT ~ Pi 
PT + Pi 



(4) 
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FIG. 1: A quantum dot (QD) connected to noncollinear fer- 
romagnetic leads, L and R, through the tunnel barriers, Vl 
and Vr, respectively. The magnetic moment in lead L(R) is 
tilted by 0/2 (—9/2) from the z axis in the z-x plane. The 
spin polarization in the leads is given by p. A single level eo 
is considered in the QD. 



(0 <p < 1). 

In the QD, a single energy level eo is considered. The 
Zeeman energy, Ez = gp^B, can be taken into account 



e o,2TM — e o ± 



E 7 



(5) 



where B is an external magnetic field in the z direction. 
The magnetic field produced by the ferromagnetic leads 
may be included in B. The Coulomb interaction for dou- 
ble occupancy in the QD, U, is assumed to be infinitely 
large. 

The strength of tunnel coupling to the leads is charac- 
terized by 



n(PT+Pi)(\v L \ 2 +\v R n 



(6) 



We define an asymmetric factor of tunnel couplings by 



(7) 



(— 1 < v < 1). In the text, we treat the case of symmetric 
barriers, Vl = Vr = V (v = 0). Then we can regard the 
z axis as a well-defined quantization axis in the QD from 
the symmetry of the system. The generalization to the 
case of Vl ^ Vr (v ^ 0) is straight-forward, as explained 
in Appendix A. 

To extract conduction modes which couple to the QD, 
we perform a unitary transformation for the conduction 
electrons in the two leads. In the energy range of [e, e + 
de], the number of majority spins, p-\de, is larger than 
that of minority spins, p\_de, by (pj — Pi)de in each lead. 
P\_de of majority spins are paired with minority spins, 



9 0) 

a feT = cos|c^ 



8 ( a ) 



(a) . . @ (s) 

a k] = cos fcj^ + sin jc\,\, 

(s) (a) 

a kl = cos |cy, - sinfcy. 
a kl = cos jc^l +sm lC y, 



(8) 



where e < e k ^, t k ' [ < e + de. The rest of majority spins 



is transformed by 



a k ] 



_ _» 
feT ■ 



a kl = -c 



(9) 



Here, we have introduced the symmetric and anti- 
symmetric combinations for the conduction electrons in 
the two leads, 



» 
(a) 



(cLfc,Z| + CRkrf ) , 
(— CLfc,2t + C_R,fe, r |) 
^ (Cifc,U + CR k , r l) , 
(cLkdl ~ CR k . r i) . 



(10) 



Note that index a =|, J, of and afe CT does not any more 
mean majority or minority spins. By this transformation, 
the tunnel Hamiltonian Ht [the last term in Eq. is 
rewritten as 



k a 

+V2V " ( cos 7 a lt^T + sm J a kl^l + n -' 



(11) 

The summations of and ^ correspond to the density 
of states pi and p-\ — pi = Ap, respectively. The modes 
a k f and a k ^ are completely decoupled from the QD and 
hence they can be disregarded in the discussion of spin 
splitting of the QD level and Kondo effect. 
The total Hamiltonian {1} is rewritten as 



H =J2 £ kaal a a ka 

k.a 

+ ^ ^,zad\ a d za + Ud\^d z ^d\^d z i 

a 

+V2Vj2( A Uz<r+h.c), 



where 



-4i 



k 

^2' a k l 



feT) 



(12) 

(13) 
(14) 



The summation over k is ^2 = + in the first term 
in Hamiltonian (fT2|) . The corresponding density of states 
is pp 

III. SCALING THEORY 

In this section, we apply the poor man's scaling 
theory^! to the Hamiltonian JT5|) and evaluate the spin 
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splitting of QD level Se and Kondo temperature. The 
scaling procedure consists of two stages^ In the first 
stage (Sec. III. A), we reduce the energy scale from the 
band width Dq to D\ where the charge fluctuation is 
quenched. By integrating out the excitations in the en- 
ergy range of D% < D < Do, we renormalize the energy 
levels in the QD, e^ and ej_, and determine the spin split- 
ting Se — ei — £•(-. Even when B = 0, we find a finite 
Se. In Sec. III. B, we proceed to the second stage scal- 
ing for D < Di where the spin fluctuation is dominant. 
We make the Kondo Hamiltonian and discuss the Kondo 
effect. 



A. Spin splitting of QD level 

We examine the Coulomb blockade region with one 
electron in the QD, — Dq <C eo *C p — A and /i + A <C 
e + U, where p, = is the Fermi energy of conduction 
electrons in the leads. The Coulomb interaction U is 
assumed to be strong enough to forbid double occupancy 
in the QD; U — > oo. 

The number of electrons in the QD fluctuates between 
zero and one through Ht- We denote |0) for the empty 
state and \za) = d\ a |0) for singly occupied state with 
spin za. They have the energies Eq and E za , respectively. 
Initially, Eq — and E Z(T = eo, ZCT . 

Under the instruction of the poor man's scaling, we 
renormalize the energies Eq and E za by integrating out 
the high energy excitations in the conduction band. Re- 
ducing half of the band width from D to D — \dD\, Eq 
and E za are renormalized as Eq + cIEq and E za + dE za , 
where 



dE Q 

dE z1 
dE zl 



2\V\ 2 \dD\ ( 2 e 

! — — ■ — pi + Apcos 2 - 

D + E z] - E V 4 

2\V\ 2 \dD\ ( . . 2 6 

2\v\ 2 \dD\ ( a 2 e . 

D + E - E zl \ Fi F 4 



(15) 
(16) 
(17) 



within the second-order perturbation with respect to Ht- 
For D ^> \E Z(7 — Eq\, we obtain the spin-dependent scaling 
equations for the renormalized energy levels, e Z(T = E za — 

Eq, 



de za _ T g 
dlnD ~ ~ tt 



where 



r T/l = 



lip COS ; 



(18) 



(19) 



is the level broadening for spin z t /z | in the QD. Note 
that the renormalization of e za depends on the coupling 
strength for the opposite spin, 



Integrating Eq. (JTHJ from (D ,e 0]Za ) to (Di,e zcr ), we 
obtain 



£za — tO.za H m — — , 

TT D\ 



(20) 



where Di sa ~^z\ a t which the perturbation theory 
breaks down^ D\ satisfies the equation 



Di « e , zT H In — — . 

7T Dl 



(21) 



The charge fluctuation is quenched at the energy scale of 
D ~ D%. As a result, the renormalized spin splitting is 
given by 

Se = e zi - e zT = — pcos^ln^ - Ez- (22) 

7T Z L>l 

This has been obtained for the P alignment (6 = 0) in 
Ref. [H. 

We observe a finite Se in the absence of magnetic field, 
Ez = 0. We estimate the value in the experimental sit- 
uation using Ceor^ where D /A — 50 and e /A = —2. 
First, Di is determined by solving Eq. (|2"Tj) . Then Eq. 
(|2"2"| with E z = yields the spin splitting 5e. Figure [H 
shows Se and Di as functions of pcos(9/2). Since D\ 
depends on pcos(9/2) only weakly, 5e is almost linearly 
proportional to p cos(0/2) (the fitting is indicated by solid 
line in Fig. 2). In the presence of spin polarization, p =/= 0, 
the spin splitting decreases with increasing 9 from to 
tt. 5e is maximal in the P alignment (9 = 0) and zero in 
the AP alignment (9 = tt). In the experimental result^ 
Se = 9meV with p — 0.31 and A = 30meV, and hence 
Se/A — 0.33 in the P alignment. The calculated result 
indicates Se/A = 0.338 when p = 0.31 and 9 = 0, which 
is in good agreement with the experimental observation. 

In the case of asymmetric barriers, Vl Vr, the spin 
splitting is calculated in Appendix A [see Eq. (|Al0j) ]. 



A / 9 6 9 • 2 . D _ 

Se = — p\ cos 2 - + v 2 sin -In— E z , 

tt V 2 z Di 



(23) 



with an asymmetric factor v in Eq. ([7])i2i In the absence 
of magnetic field (Ez = 0), Se is always maximal at 9 = 
and minimal at 9 — tt although Se ^ in the AP 
alignment (9 = tt) in general. The observation of a finite 
spin-splitting with AP alignment— may be explained by 
Eq. (|23|) with asymmetric barriers (v ^ 0). 



B. Kondo effect 

In the energy scale of D <C D\, the empty state |0) is 
irrelevant: The number of electrons in the QD is fixed 
to be (rid) — 1- To restrict the QD states to \z f) or 
\z I), we apply the Schrieffer- Wolff transformation to the 
Anderson model with renormalized levels, e 2 | and e z ±, 
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=5 



0.4 0.6 

pcos(0/2) 



FIG. 2: Spin splitting Se = e z | — e z -\ of the QD level (open 
circles) as a function of pcos(0/2), when Do/ A — 50, eo/A = 
—2, and Z?z = 0. The tunnel barriers are symmetric (v = 0). 
<5e is fitted to a straight line, Se/ A = 1.046p cos((9/2) (solid 
line). Closed circles indicate Di which is determined by Eq. 
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When Z) <C <5e, state |z J.) is irrelevant and hence the 
spin-flip processes do not take place. The coupling con- 
stants cease to increase with decreasing D. 

Let us consider a situation where Se is negligibly small 
(Se -C Tk). This situation can be realized by tuning 
an external magnetic fields Then the scaling equations, 
([27)) and ([28]) . are applicable until the scaling reaches 
a fixed point of strong coupling limit. As a result, the 
Kondo effect takes place. The initial condition for the 
scaling is given by 



J± = J; 



2? 



J. 



2V 2 



= J 



(29) 



with D = D\ 1 where e = e z \ = e z [. e = —D± from the 
definition of D\ in the previous subsection. 

From Eqs. (f27|) and |28|) . we find the scaling trajectory, 
(k^+ktf-Akl = (p 1 +p l ) 2 J 2 p 2 cos 2 (9/2) and fc T -fc; = 
G°T + Pl)Jp C0S (®/2)- O n the trajectory, the coupling 
constants go to the fixed point of k±, fc-f, fc[ = oo, with 
decreasing the energy scale D. The energy scale where 
the fixed point is reached yields the Kondo temperature, 



given by Eq. (|2"0|) . We obtain the Kondo Hamiltonian 



#Kondo =Hq 



k.a 



sd 



--J + S+A\At + J_S~A\A i 
+ S z (J z] A\Ai-J zl A\A l ), 



(24) 
(25) 

(26) 



where S's are spin operators in the QD, S + = d z +d z i, 

S- = d} zi d zh and S z = (4 T d»t ~ 4l d zl)/ 2 - Here ' 
we have neglected a small level shift, — £q[(p-f + 
J2<j Jza\^ + &P cos(8 /2)\d\ a d za , and potential scat- 
tering terms: The former may be included in the second 
term in Hq, whereas the latter is not relevant to the 
Kondo effect. J+ = J_ = J± and J za are exchange 
coupling constants. They are given by 



J± = V 2 



J z 



2V 2 



We continue the scaling procedure to D < Di, now 
using the Kondo Hamiltonian (|2"4"|) . We obtain the scaling 
equations for Jj_ and J za to the second-order with respect 
to 7J s( j when D 3> Se. 



dk± 
dlnD 

dk ° _ 91,2 



(27) 
(28) 



where k± and fc<j are effective coupling constants, k± 
[(PT + Pi )/2]J ± Vl-P 2 cos2(0/2) and = [(p T 
Pi )/2]J ZCT [1 + op cos(0/2)]. 



Tk(Q,p) = Tqexp 



arctanh (pcos |) 



pCOS 2 



(30) 



This is an extension of the result in Ref. |l3| for the P 
alignment i 2 ^ 

Equation ([30)1 describes the dependence of the Kondo 
temperature on the spin polarization p and relative an- 
gle 9 between the magnetic moments in the ferromagnetic 
leads. Tk(0,p) is a function of pcos(9/2). With an in- 
crease inpcos(0/2), Tk decreases as shown in Fig. 3 (solid 
line), reflecting the suppression of spin fluctuation in the 
QD. When the leads are completely polarized (p = 1) 
and in the P alignment (9 = 0), the Kondo effect vanishes 
owing to the absence of spin fluctuation. In the AP align- 
ment (9 = 7t), Tk is independent of the spin polarization 
p and given by Tk = D\ exp{— l/[(pj + pi)J}}. This 
expression coincides with that for non-magnetic leads 

( P = 0)M 

The properties of this Kondo effect are determined by 
the fixed point of strong coupling limit. Although J±, 
J z |, and J z i —* oo for D — > Tk, their ratios are finite, 



Jz\ 1— PC0S : 
Jzi 1+pCOSi 



J 



* J. 



The former expresses an important characteristic of the 
spin asymmetry by the coupling to noncollinearly spin- 
polarized leads. The ratio of J z \/J z i is equal to that of 
effective density of states of conduction electrons coupled 
to \z I) and \z |), pj/pf, where p\ = p l + Aps\n 2 (9/A) 
and p* = Pl + Apcos 2 (6>/4) [see Eqs. and ([T4])]. In 
other words, the relation of p*yJ z i — P* z [Jzi is realized 
in the strong coupling limit. In the formation of Kondo 
singlet state, spin state \z t ) in the QD is coupled to 
conduction electrons A± with density of states p*^ by J z i, 
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the constraint 



0.4 0.6 

pcos(6/l) 



FIG. 3: The Kondo temperature Ik obtained by the scal- 
ing calculation (solid line) and slave-boson mean-field theory 
(dashed line), as a function of pcos(#/2). The tunnel barriers 
are symmetric (v = 0). The spin splitting is Se = 0. Tk is 
normalized by the value at pcos(#/2) = 0. In the scaling cal- 
culation, we choose (pf + Pi)J = A/(7r|e|) with renormalized 
QD level e = -Di = -1.43A3 



whereas spin state \ z J.) is coupled to conduction electrons 
Af with density of states p* z ^ by J z |. The relation of 
p* z ^J z -\ = p*|J z | indicates that the spin-up and -down 
in the QD are equivalently screened and a local singlet 
state is realized in spite of the spin polarization in the 
ferromagnetic leads. 

In the case of asymmetric barriers, Vl ^ Vr, the 
Kondo temperature Tx(0 7 p) is given by Eq. (|A11[) in 
Appendix A. The above-mentioned discussion is appli- 
cable to the case of v ^ if pcos(#/2) is replaced by 



pycos 2 ((9/2) +v 2 sin 2 (6»/2). 

We have observed the Kondo effect for Se <C Tk. In the 
opposite situation of Se S> T K , the scaling by Eqs. (j2"T|) 
and (f2"5| ends at D ~ Se and then the coupling constants 
saturate. No Kondo effect is expected in this case. 



IV. SLAVE-BOSON MEAN-FIELD THEORY 

From the scaling calculation, we obtain an analytic ex- 
pression of the Kondo temperature, Eq. (|30p. when the 
spin splitting Se <C Tk- The Kondo effect does not take 
place when Se 3> Tk- To evaluate the Kondo tempera- 
ture with Se ~ Tk, we adopt the slave-boson mean-field 
(SBMF) theory^ in this section, assuming that U = oo. 
The conductance is also calculated at temperature T = 0. 

In the slave-boson formalism, we introduce an auxiliary 
boson field to represent an empty state in the QD. Using 
the creation operator of the boson b\ the empty state 
in the QD is written as |0) = 6^|vac), where |vac) is the 
vacuum state. The singly-occupied states are expressed 
as \za) = /] CT |vac) with a =t,4, where j\ a is a pseudo- 
fermion operator to create the spin state. When U = oo, 



fzafza 1 



(31) 



is required since the QD state should be empty or singly- 
occupied with either spin. The original operators for an 
electron in the QD are rewritten as d\ a = fj a b and d za — 

We apply this formalism to the Anderson model ([TJ) 
with renormalized QD levels, Eq. (|2H)) . obtained by the 
first stage scaling in Sec. III. A. Hence the charge fluc- 
tuation in the energy range of D\ < D < Dq has been 
taken into account. The density of states is constant, 
Pyi, in the band of — D± < u> < D±. The Hamiltonian is 
rewritten as 

H = E e k*C LkM C Lk ,lc + ^2 e k<7 C Rk ra C Rk ,ra 



ka 



ha 



+ Y,( Vc UlJ^ + Vc Rk,rJr^ + h.C.), (32) 



where A is a Lagrange multiplier to consider the con- 
straint, Eq. (f3"Tj) . In the following discussion, we fix D\ 
and e = (e zA + e zl )/2. 

In the SBMF theory, the Kondo effect takes place when 
the boson field is condensed: The boson operators b, 
are replaced by a c number, (b) = r. Then the Hamilto- 
nian is reduced to a non-interacting Anderson model in 
which the QD levels are given by E z ^ — e z ^ + A and 

the tunnel coupling is V — rV ~ We can easily calculate 
the Green functions in the QD 



G zazcr (to) 



1 



■iT„ 



G 



= G zLz] (uj) = 0, 



(33) 
(34) 



where T a = r 2 T a . The fact that the z axis is a well- 
defined quantization axis results in zero off-diagonal ele- 
ments. Minimizing the expectation value of the Hamil- 
tonian with respect to A and r, self-consistent equations 
for them are obtained as 



7T £>l 

tT 

^ + (n f ) = 1, 



where (nf) — (n z ^) + (n z i), with 



= {ftofzo) = - arctan 



7T 



(35) 
(36) 

(37) 



We are interested in the case of €20- <C p — A, where 
(nf) ~ 1 (Kondo regime). By solving Eqs. (f3l>|) and (pirjj) . 
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FIG. 4: The Kondo temperature 2k obtained by the slave- 
boson mean-field theory as a function of the spin splitting 
Se. The tunnel barriers are symmetric (y = 0). The spin 
polarization in the ferromagnetic leads is (a) p = 0.1, (b) 0.5, 
(c) 0.9, and (d) 1. The angle 9 between the magnetic moments 
in the leads is 0, 7r/2, and n. Both 7k and Se are normalized 
by T K (Se = 0,6 = 7T,p). Note that 7k = for any Se in the 
case of (d) p = 1 and 9 = 0. 



we observe the spin-split Kondo resonant levels, which 
are located at 



Se , 

E *Vzl = =Fy I l±pcos- 



(38) 



below and above the Fermi level \i = 0. Note that they 
are not equi-distant from \x = except in the AP align- 
ment [9 = 7r). The resonant widths are given by r-m. 
We define the Kondo temperature by the geometric mean 

of the resonant widths, Xk = 



r T ri. 



It is given by 



T K (Se,9,p) = JT K (0,6,p) 2 - \E z] E. 



(39) 



and 



T K (O,0,p) = £>iexp 



— : V cos — arctanh v cos — 

A 2 V 2 



(40) 



In Eq. gDJ, Se = and thus e = e A = e zl (= -D x ). 

First, we compare Eq. (|40|) with Eq. (|30|) obtained by 
the scaling theory when Se = 0. In the case of non- 
magnetic leads (p — 0) or of magnetic leads in the AP 
alignment (9 — n), they are identical to each other: 
T K = Diexp(-7r|e|/A) = D x exp{-l/[(p T + Pi ) J]}. In 
the case of magnetic leads, both of them are functions 
of pcos(0/2). They decrease with increasing pcos(9/2) 
similarly, as shown in Fig. [3] The Kondo effect is 
the strongest at pcos(9/2) = 0, whereas it disappears 
pcos(9/2) — 1. This semi-quantitative agreement clearly 
indicates an applicability of the SBMF theory to the 
present problem. 



Second, we discuss the Kondo effect with finite Se. In 
Fig.Hl we plot T K (Se, 9,p) in Eq. (J35J) as a function of Se. 
It decreases monotonically with an increase in Se, reflect- 
ing the separation between the Kondo resonant levels for 
z | and z I in Eq. (|38[) . The Kondo temperature van- 
ishes at Se = 2T K (0,9,p)/y/l-p 2 cos 2 (6»/2). Note that 
this sudden disappearance of the Kondo effect is an ar- 
tifact by the SBMF theory. The theory works unless the 
Kondo temperature is too small^ 

The suppression of the Kondo effect can be also un- 
derstood in terms of the spin accumulation, (Snf) — 
(n z f) — (n z i). Using Eq. |37|) and (rif) ~ 1, the Kondo 
temperature is rewritten as 

T K (Se) ~ T K (0)^Jsm (n(n z1 )) sin (ir(n zi )) 



T K (0)cos ( -(Sn f ) 



(41) 



As the spin accumulation (Snf) increases, the spin fluc- 
tuation is suppressed, which weakens the Kondo effect. 
T K {Se) ~ when (Sn f ) ~ 1 {(n z] ) ~ 1 and (n zL ) ~ 0). 
In this situation, the Kondo effect is very weak because 
spin z I hardly exists in the QD. 

The above-mentioned discussion can be generalized to 
the case of asymmetric barriers, Vl Vr. The Kondo 
temperature Tk(0, 9,p) is given by Eq. (|A12|) in Ap- 
pendix A. It is obtained by replacing pcos(9/2) in Eq. 



PTJf to be pycos 2 (6>/2) + v 2 sin 2 (9/2). 

Now we discuss the conductance G through the QD in 
the case of symmetric barriers. The SBMF calculation 
yields 



G = 



2(1 -P 2 ) 
1 — p 2 cos 2 

T K (5e,9,p) 



■ p 2 sin 2 



Se 



2 \T K (O,0,p) 



T K (0,9,p)J ■ (42) 

The conductance is explained as the resonant tunnel- 
ing through the Kondo resonant levels at E z y z ^ in Eq. 
(|38f with the width of r-m. With non-magnetic leads 
(p = 0), the conductance is simply given by G = 
(2e 2 /h)[l - {Se/2T^ } ) 2 }, where is the Kondo tem- 
perature at Se = 0. 31 When Se = 0, the resonant tun- 
neling through the level of E z y z ^ — yields the unitary 
limit, G = 2e 2 /h. A finite Se splits the resonant level 
into two, which reduces the conductance. At Se > 2T-^\ 
G = since the Kondo effect does not take place in the 
SBMF theory. With ferromagnetic leads (p ^ 0), we 
expect a suppression of G below G = 2e 2 /h owing to 
the spin- valve effect j§j§iL2aL besides the splitting of the 
Kondo resonant levels by finite Se. 

First, we consider the case of Se = to discuss the 
spin-valve effect only. Then the Kondo resonant levels 
are fixed at the Fermi level, E z y z ^ = 0. Equation 
reduces to 



G = 



2e 2 



l-p 2 



h 1 



p z cos^ 



(43) 
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FIG. 5: Conductance G through the QD as a function of 6, in 
the absence of spin splitting (Se = 0). The spin polarization 
p in the leads is varied from to 1. The tunnel barriers are 
symmetric (v = 0). 



G is plotted as a function of 9 in Fig. [5J for various p. 
When < p < 1, the conductance decreases with an 
increase in 9. This is due to the spin- valve effect in which 
the conductance is maximal in the P alignment (9 = 0) 
and minimal in the AP alignment (9 — n). Change of 
G, G(9 — 0) — G(9 = 7r), is more prominent for larger p. 
When p = 1 (half-metallic leads) , the conductance always 
vanishes for the following reason. At 9 = 0, this is due to 
the absence of Kondo effect. At 9 ^ 0, this is ascribable 
to an interference effect between two spin-components. 
This is analogous to the Fano-type antiresonance in a 
quantum wire with a side-coupled quantum dot i 32 ' 33 i 34 
In Appendix B, we examine a non-interacting model in 
which a quantum dot is coupled to two ferromagnetic 
leads with perfect polarization (p = 1). We find a dip of 
G (G — 0) when an energy level in the QD matches the 
Fermi level in the leads. The dip can be understood as 
the Fano-type antiresonance because one spin component 
in the QD is coupled to both the leads (corresponding 
to "a quantum wire") while the other spin component 
is coupled to one of the leads ("side-coupled quantum 
dot" ) . In the case of p = 1 in Fig. [5J the dip makes 
G = since the Kondo resonant level always appears at 
the Fermi level. 

Second, we examine the conductance G in the presence 
of spin splitting Se. In Fig. [51 G is plotted as a function of 
Se, for various p and 9. Se in the abscissa is normalized by 
Tk(5e = 0, 6 = 7r,p), the Kondo temperature at Se = 
and 9 = it for each value of spin polarization p in the 
leads. When p is small [(a) p = 0.1, (b) 0.5], G decreases 
with increasing Se and becomes zero at a critical value 
of Se. This is due to the splitting of the Kondo resonant 
levels for z f and z |, as in the case of non- magnetic 
leads. When p is large [(c) p — 0.9, (d) 1] and 9^0, the 
second term in Eq. (|42|) makes a non-monotonic behav- 
ior of G(5e). At p ~ 1, the conductance is suppressed by 
the dip of Fano-type antiresonance at Se = 0. With in- 
creasing Se, the dip is shifted from the Fermi level, which 
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FIG. 6: Conductance G through the QD as a function of spin 
splitting Se. The spin polarization in the leads is (a) p — 0.1, 
(b) 0.5, (c) 0.9 and (d) 1. The angle 8 between the magnetic 
moments in the leads is 0, n/2, and ir. The tunnel barriers are 
symmetric (v = 0). Se is normalized by Tk(5e = 0, = 7T,p). 
Note that G = for any Se in the case of (d) p = 1 and 9 — 0. 



increases the conductance. This results in a maximum of 
G as a function of Se, in combination with the suppres- 
sion of G by the splitting of Kondo resonant levels. 



V. CONCLUSIONS 

We have theoretically investigated the Kondo effect in 
a QD coupled to noncollinear ferromagnetic leads, using 
the poor man's scaling method and the SBMF theory. 

We have used the poor man's scaling procedure in 
two stages. In the first stage, we have reduced the en- 
ergy scale D until the charge fluctuation is quenched 
and renormalized the energy levels for spin | and J, in 
the QD. We have found a finite spin splitting Se — 
£j — e-f in the absence of magnetic field. This stems 
from the spin-dependent charge fluctuation. The spin 
splitting takes place in an intermediate direction be- 
tween the magnetic moments in the two ferromagnets. 



Se cx p\j cos 2 (#/2) + v 2 sin 2 (0/2), where p is the spin po- 
larization in the leads, 9 is the angle between the mag- 
netic moments, and v is an asymmetric factor of tunnel 
barriers (—1 < v < 1). Hence the spin splitting is al- 
ways maximal in the parallel alignment of two ferromag- 
nets (9 = 0) and minimal in the antiparallel alignment 
(9 = it). 

We have proceeded to the second stage of the scal- 
ing for D < D\, where a localized spin fluctuates in the 
QD with fixed spin splitting Se. When Se is irrelevant 
(Se <^ Tk), we have derived an analytical expression of 
the Kondo temperature Xk(#,p). Tk is a decreasing func- 
tion with respect to p*J cos 2 (#/2) + v 2 sin 2 (#/2). When 
Se 3> Tk, the Kondo effect does not take place. 



We have applied the SBMF theory to the renormal- 
ized Hamiltonian by the first stage scaling in which the 
charge fluctuation in the energy range of D\ < D < D 
(original band width) has been taken into account. When 
Se = 0, the Kondo temperature decreases with increasing 

pj cos 2 (#/2) + v 2 sin 2 (#/2), in semi-quantitative agree- 
ment with that obtained by the scaling method. This 
clearly indicates an applicability of the SBMF theory to 
the present problem. When Se is comparable with Xk, we 
have observed a splitting of the Kondo resonance, sup- 
pression of the Kondo effect, and spin accumulation in 
the QD. Tk{S€,6,p) has been evaluated by the width of 
the split Kondo resonances. We have also calculated the 
conductance G at temperature T = 0. G could show 
a non-monotonic behavior with increasing Se, in contrast 
to 7k. This is ascribable to an interference effect between 
two spin components, analogous to the Fano resonance. 
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APPENDIX A: CASE OF ASYMMETRIC 
BARRIERS 

In the case of asymmetric barriers, Vl ^ Vr, we can 
choose the quantization axis for the QD spin, as in the 
following. We define an axis n which is tilted by ip from 
the z axis in the z-x plane as shown in Fig. [7] We denote 
n | /to I for spin- up/down in the direction. Fermion 
operators in the QD, di a , d ra , are related to d na by 




FIG. 7: Quantization axis n for the QD spin in the case of 
asymmetric barriers, Vl ^ Vr. l(r) indicates the direction of 
magnetic moment in lead L(R), which is tilted by 9/2 (—9/2) 
from the z axis in the z-x plane (see Fig. 1). The axis n is 
tilted by ip from the z axis in the same plane. 



The coefficients are given by 



r m 

A 2 
A 3 



cos I cos % + v sin | sin ■ 
VL/jsinf sinf , 



(A5) 



^ A 4 = v LR sin | 



Cka and 5^ are defined as 



4 cosf. 



c k] — Tj(^L c Lfc,iT + VR c Rk,rl) 
— I ( M_ „_ . . . i T/L „_ _ 



Cfc-[ 
CfcJ. 
CfcJ. 



VRCRk,r l)t 
V~LCRk,rl), 



(A6) 



where V = y/V£ + VJ, v LR = 2V L Vr/V 2 , and v is the 
asymmetric factor defined by Eq. l[7|). Then the tunnel 
Hamiltonian Ht [the last term in Eq. (JTJ] is rewritten as 



(A7) 



where 



dii 



% -^d 



a \ 



cos 



^d 

V-2ip i 
4 a nJ j 



d r1 = cos ^j^d n1 - sin ^±pd n |, 



e -+kd 



cos e -±^d nl . 



(Al) 
(A2) 



We consider a following transformation for conduction 
electrons with e < efcf, £fc'^ < e + de. For /9jde of majority 
spins, 



afeT = Aic feT - A 2 c feT + A a c k/l - A4CV1, 
0*4 = -A 3 c*;| - A 4 c feT + AiCfc/j + A 2 c*/;, 
5*4 = A 2 c k -\ + Aic feT + A 4 c fe ^ + A 3 c k 'i, 
0*4 = AiCk-[ - A 3 Cki - A 2 c k >i + A1CV4. 

For the rest of the majority spins, 



(A3) 



.(AiCfcf - AgCfcf), 

■■(-A 3 Cki - A 4 c feT ), 



(A4) 



^1 = Efc 'flfel 



v/aT 
VaT 



_A|E fe "^T, 

AtXVW 



(A8) 



The modes of 6*4 and 0*4 are decoupled from the QD. 

Although Eq. (|A3[) is the unitary transformation with 
arbitrary coefficients, Eq. (|A4p is not. From the condition 
[A|,Aj] + = 0, that is, AiA 3 = A 2 A 4 , we determine ip 
(direction of axis n) so that Eq. (|A4[) becomes the unitary 
transformation. It is given by 



tan if = v tan ■ 



(A9) 



The quantization axis n leans toward the direction of 
magnetic moment with larger tunnel coupling from the z 
direction. 

We calculate the spin splitting, Se = e n ^ — e„t, in 
the same way as in Sec. III. A. In the most right- 
hand side of Eq. (|22[) . pcos(9/2) should be replaced 
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by p cos(6> / '2)/ cos <p = pJ 'cos 2 (0/2) + v 2 sin 2 (<9/2). This 
yields 



<5e = -m/cos 2 ^ + 1; 2 sin 2 £ In ^ - E z . (A10) 
7r V 2 2 D\ 

In the case of asymmetric tunnel barriers (v ^ 0), we 
observe a finite spin splitting even in the AP alignment, 
5e = (A/tt)p|u| ]n(D /Di), when E z = 0. 

The Kondo temperatures are obtained using the scal- 
ing method (Sec. III. B) or by the SBMF theory (Sec. 
IV). Replacing pcos(6>/2) in Eqs. ([3D]) and @J) by 



p\l cos 2 {9/2) + v 2 sin (9/2), they are generalized to 

1 



T K (6,p) = J Diexp 



arctanh ( p*/ cos 2 § + v 2 sin 2 | 



P\l cos 2 | + v 2 sin 2 | 



(All) 



and 

T K (O,0,p)=£) 1 exp 



A 



-p^ cos 2 - + v 2 sin 2 -arctanh ^py cos 2 - + v 2 sin 2 - 

(A12)~ 

respectively. With an increase in the asymmetric fac- 
tor I?; |, the spin fluctuation in the QD is suppressed and 
accordingly the Kondo effect is weakened. 



APPENDIX B: NON-INTERACTING MODEL 
WITH p = 1 

In Appendix B, we examine a non- interacting model 
in which a QD is coupled to noncollincar ferromagnetic 
leads. We assume that the leads are half-metals (pi = 0) 
and thus the spin polarization is p = 1. With U = 0, the 
Hamiltonian (TTJ reads 

H = ^ efc T c L/MT CLfe >'T + ^2 ek l c Rk,r\ c Rk,ri 
k fc 

+ ^e 4r<4 CT + H T . (Bl) 

cr 

Here, the tunnel Hamiltonian is written as 

H T = Y.^UlA* + V * C Rk,rJr* + h.C.), (B2) 



or 



H T =Vj2i A Una+h.C.), 



(B3) 





1 


- 6=0 — , 










e=m — - / 




(a) 






6=71/1 //'\ ;' 








6=37t/4 // j j 


v\ 








6=7i —■-// \ \ 








0.5 


// ■■■■ '• 1 




_ 


o 




yy /' \ ii 

^5*' •**' \ !! .' 








n 


•***" '■ '■! ■' 
' '■■■):,■ 

- - - : - - - 








In 




\ 












(b) 




3n/2 










\ \ 1 
V \ 






© 


% 








eo 






y 






■K/l 




'\'% 









1 1 


i 





-2-10 1 2 

eo/A 

FIG. 8: Calculated results using a non-interacting model in 
which a QD is coupled to two noncollinear ferromagnetic leads 
with perfect polarization (p = 1). (a) Conductance G and (b) 
phase shift at the Fermi level <5?7(0), as functions of QD level 
eo- eo is normalized by the tunnel coupling strength A given 
by Eq. The tunnel barriers are symmetric (v = 0). The 
angle 6 between the magnetic moments in the leads is shown 
in panel (a). 



using the quantum axis n in Eq. (|A9|) in Appendix A. 

By virtue of the quantum axis n, the Green functions 
in the QD can be written in a diagonal form, 



Gr, 



i 



lu - e + iT na 

Gn]nl{u) = G n |„f (cj) = 0, 



(B4) 
(B5) 



where 



nl/nl 



A 
~2 
A 
~2 



1 ± 



cos(0/2) 

COS If 



l±^Jcos 2 (9/2)+v 2 sin 2 (6>/2) 



(B6) 



are the level broadenings for spin n ] /n { in the QD. 
The T-matrix T through the QD is given by 



{Rk,r]\f\Lk,l])=V L V R 



cos ■ 



2<p 9 + 2<p ^ 

cos : — G nT „ T (u;) 



sm ■ 



4 4 
9- 2(p . 9 + 2ip 



4 4 
VlVr(lj - e ) cos I 



Gnlnl(u) 



(uj - e + ir„|)(w - e + ir„|) 



(B7) 
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with to being the energy of incident electrons. Using the 
T-matrix, we obtain the conductance G through the QD 



G= -(2tt Pt 



2 \(Rk,r^\f\Lk',n)\ 2 
e 2 4r L r fl e 2 cos 2 f 



LO=0 



(B8) 



= np^V 2 and Vr = np^V^ are the strength of the 
tunnel coupling to leads L and R, respectively. A = 
T l +Tr in Eq. © for Pi = 0. 

In Fig. EJa) , we present the conductance G as a func- 
tion of QD level eo. In the P alignment (6 = 0), r„T = A 
and r ni = in Eq. ([56]) . Then Eq. JB8} is reduced 
to G = (e 2 /h) ■ 4r L r it /(eg + A 2 ), which indicates the 
usual resonant tunneling through QD level eo by one spin- 
component [solid line in Fig. [5[a)]. When < 9 < tt, 
however, Eq. (|B8[) yields G = at eo = 0. We observe 
a dip of G at eo = in Fig. E^a). I n the AP alignment 
(9 = 7r), G = for any value of eo since we are considering 
half- metallic leads with pi = 0. 

The dip of the conductance at < 9 < tt is attributable 
to an interference effect between two spin components. 
Consider the tunnel coupling between lead L and QD in 
tunnel Hamiltonian (|B2|1 . The QD state \l f) is coupled 
to the lead L by Vl, whereas |^ J.) is decoupled since p[ = 
in the lead. On the other hand, both the QD states |Z f) 
and \l [) are coupled to the lead R by Vrcos(0/2) and 
— Vrsui(#/2), respectively, because the QD state \r |) is 
coupled to the lead R by Vr and the QD state is rewritten 
as |r t) = cos(6>/2)|Z f) - sin(0/2)|J |). This situation is 
analogous to that of a quantum wire (lead L - QD state 
|/ |) -lead R) with a side-coupled QD (QD state \l 1) 
-lead i?) , in which the Fano-type antiresonance has been 
observed . 33 ' 34 In general, an interference between a local- 
ized state and a continuum of states results in an asym- 
metric Fano resonance, G oc (io + q) 2 / (e§ + l), where eo is 
the resonant level normalized by the resonant width. 32 In 
a geometry of a quantum wire with a side-coupled QD, 
the Fano factor is q = which makes a dip at eo = 0. 

An alternative explanation for G = at eo = is given 
in terms of phase shift. Using the T-matrix in Eq. (|B7[) . 



the phase shift 5rj(uj) is defined by 

5t](lu) = axg{Rk,r1\f\Lk,l-\) 

(ep -cj)(r» T + r ni ) 
(e - uj) 2 - r„ T r ni ' 



arctan 



(B9) 



Using Eq. (|B6|) . the conductance G in Eq. (|B8[) is rewrit- 
ten as 



G 



e 2 4r L r fl 

h{T L +T R ) 2 " 2 



cos 2 -sin 2 ^(O), (BIO) 



where 8rj(0) is the phase shift for electrons at the Fermi 
level, u) = 0. 

Next, we express the number of electrons in the QD 
(rid) in terms of the phase shift Sr/(0) (Friedel sum rule). 
Using the local density of states in the QD 



Pd 



we obtain 



(uj) = - V* -lmG ncmrT (uj) 

1 ' TT 



1 d8ri(oj) 

TT du> 



(rid) = / Pd(u)duj = . 



(Bll) 



(B12) 



This sum rule is different by a factor of one half from the 
conventional Friedel sum rule in the case of non-magnetic 
leads. 

Figure b) shows the phase shift at the Fermi level 
Sr](0), as a function of QD level eo. At eo = 0, Sr](0) = 
ir/2 when 9 = and 5r](0) = n when 9^0, corre- 
sponding toG= (2e 2 /^)-4r L r i? /(r L + r i? ) 2 andG = 
in Fig. [5Ja), respectively. This can be explained using 
the sum rule of Eq. (|B12p . In the half-filling case of 
eo = 0, the number of electrons should be unity when 
9^0. Then Eq_JHj2} yields Sr](0) = tt and hence 
G = by Eq. (|B10j) . In a special case of 9 = 0, 
the QD state \z J.) is always empty since it is com- 
pletely decoupled from two half-metallic leads. Then 
(rid) = 1/2 at eo = and thus 5i](0) = ir/2, which leads 

toG = (e 2 // l )-4r L r fl /(r L + r fl ) 2 . 
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